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Evidence for Cold Accretion onto a Massive Galaxy at 
High Redshift? 
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ABSTRACT 

In this letter we report on the discovery of a 2=2.83 Lyman-a Blob (LAB) found 
in our wide field narrow-band survey within the Spitzer First Look Survey region. The 
blob is extended over at least 95kpc and has a total Lyman-a luminosity of 2.1 xlO 44 
erg s _1 . It is only the sixth LAB known of this scale (> 50 kpc), and is associated with 
an embedded continuum source in g', R, i', K & 4.5 nm bands. The LAB's optical 
spectrum shows clumpy structures and tantalising hints of a sharp red cut-off and 
shear within the Ly-a emission line. Studies of the LAB's surface brightness profile 
and of the continuum counterpart's spectral energy distribution (SED) indicate that 
the profuse Lyman-a emission is consistent with being powered by cold gas accreting 
onto a massive dark matter halo. 
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1 INTRODUCTION 

Recently, much success has been gained in picking out high 
redshift galaxies using a narrow-band filter to scan their Ly- 
man a emission in a small redshift "slice" of the distant 
universe. Steidel et al. (2000 - hereafter S00), for example, 
discovered a 'spike' in the redshift distribution of a sample 
of LBGs at z=3.09 and observed the region with a narrow- 
band filter tuned to pick out Lyman a emitting galaxies at 
that redshift to confirm the presence of a proto-cluster. They 
not only found evidence of a proto-cluster, but also discov- 
ered two very extended (~ lOOkpc) and extremely luminous 
{Lhya ~ 10 44 erg s _1 ) regions of Lyman a emission with 
very faint continuum and no radio detections. They named 
these objects "Lyman a Blobs" (LABs) on account of their 
almost monochromatic emission and amoebic appearance. 

Perhaps the most mysterious aspect of LABs is the ap- 
parent lack of an obvious source of ionisation to power such 
profuse Lyman a emission. A variety of models have been 
suggested to explain this apparent omission; three of the 
most plausible explanations are: (i) the Lyman-a emission 
comes from a dust-enshrouded, extreme starburst galaxy 
with a large scale superwind (e.g. Taniguchi & Shioya 2000; 
Ohyama et al. 2003; Wilman et al. 2005), (ii) LABs con- 
tain hidden QSOs (e.g. Haiman & Rees 2001; Weidinger, 
M0ller & Fynbo 2004; Weidinger et al. 2005) or (iii) that 
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we are observing the cooling radiation of a collapsing proto- 
galaxy inside a dark matter halo's gravitational potential 
(the so-called "cold accretion" model - e.g. Haiman, Spaans 
& Quataert, 2000, Fardal et al, 2001, Matsuda et al, 2004, 
Nilsson et al., 2006) . 

Currently only ~ 5 of the most extended (>50kpc scale) 
LABs are known (see S00, Francis et al., 2001, Dey et al., 
2005, Nilsson et al., 2006), and it is clear that more are 
required before we can pin down the source of ionization 
with any certainty, or explain the clustering or evolutionary 
aspects of LABs within the ever-increasing compendium of 
high redshift objects. 

In the following sections, we present our wide field sur- 
vey, give an account of the discovery of this new LAB, dis- 
cuss our results & state the conclusions we have drawn. 
Throughout this letter the AB magnitude system is used 
(Oke & Gunn, 1983), and a standard cosmology is assumed 
in which H = 71 km s~\ fi M = 0.27 and fi A = 0.73. 



2 OUR SURVEY 

In an attempt to pick out a statistically meaningful sample 
of z ~ 3 LABs, we have adopted a different approach from 
other narrow-band surveys which have tended to concentrate 
on covering a small or moderate sky area to large depth (e.g. 
Matsuda et al., 2004, Nilsson et al, 2006 and Saito et al., 
2006). By exploiting the comparatively large field of view 
of the Wide Field Camera on the Isaac Newton Telescope 
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(INT-WFC) our survey has concentrated on covering a very 
wide area (~ 15 sq. degrees) to relatively modest flux limits. 
In this way we can compensate for the necessarily small co- 
moving volume probed by a narrow-band filter and as a 
result we are sensitive to the brightest, most extended and 
rarest LABs. 

We have observed three areas of sky which allow us 
to take advantage of high quality public survey data. The 
fields we have observed are; the SWIRE region within the 
XMM Newton Large Scale Structure survey field (XMM- 
LSS) centred on 02 h 18 -05°00 (see Lonsdale et al., 2003); the 
Lockman Hole (LH; 10^50 +57° 00; see Lonsdale et al., 2003) 
and the Spitzer First Look Survey (FLS; 17 h 18 +59°30, see 
Lacy et al., 2005) 

Our survey comprises observations from three different 
narrow-band filters, along with complementary observations 
over the same fields in Sloan g' and i' bands. The narrow- 
band filters used (Hell, H/3b, and [OIII]r) come from the 
WFC filter set, and are centred on 4686, 4861 & 5008 A, with 
full-width half maxima of 100, 170, & 100 A, respectively. 
For the narrow-band observations discussed in this letter, we 
have used the Hell narrow-band filter; the data are found to 
be complete to HellAB ~ 24.6, g' as ~ 25.0, and Vab ~ 24.7 
respectively (5ct), while the seeing in the narrow-, g'-, and 
i'-band data are 1.03", 1.25" and 1.17" respectively. 

Our data were de-biassed, trimmed, flat-fielded, bad 
pixel- and gain-corrected, astrometrized, stacked and pho- 
tometrically calibrated using the Cambridge Astronomical 
Survey Unit WFC pipeline (Irwin & Lewis, 2001). Gen- 
eral catalogues were generated from matched pointings using 
SExtractor (Bertin & Arnouts, 1996) in double image mode. 
Sources of more than 5 adjacent pixels above la from the 
measured mean in a 4" aperture centred on the peak pixel 
are source extracted. LAB candidates are selected from our 
matched SExtracted g' and Hell catalogues based on their 
Equivalent Widths (EWs), which are estimated from their 
photometry by defining an R parameter, where 

R = io°' 4x(g AB- HcII AB) i n\ 

(from Venemans et al., 2003). Here g'AB and HellAB repre- 
sent the broad- and narrow-band magnitudes of a source. We 
can then estimate the observed equivalent widths (EW ob s) 
of each source using 

EW ob3 = 1 R ~ 1 R , (2) 

AX nb AXg 

where AA n & and AA 9 are the widths of the narrow- and 
broad-band filters respectively. LAB candidates are then se- 
lected if they have an EW obs greater than 200A and 24.1 > 
HellAB > 19.5. It should also be noted that Hayes & Ostlin 
(2006) suggest that equivalent widths estimated in this man- 
ner may be a factor of ~2 greater than the true value, how- 
ever our follow-up observations do not show this discrepancy 
(see section [3}. 

Our narrow-band data are also sensitive to low-redshift 
line emitting galaxies. The main contaminant population 
will be [OII]3727 at z ~ 0.26. Thus, with our selection of 
line emitting galaxies at EW b s > 200A, any [Oil] emission 
line would have to have a rest-frame equivalent width of > 
160 A. Hogg et al.(1998) surveyed 375 [Oil] emitters within 
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Table 1. Photometric data for the LAB. Magnitudes are mea- 
sured in apertures with 4" radii, and where they are limits, they 
are to la. The large error in the Hell narrow-band filter is due to 
several effects, the faintness of the LAB, the difficulty in fitting 
the sky background due to the very extended nature of the LAB, 
and the presence of the nearby low-z interloper to the North and 
East of the LAB counterpart galaxy. 

a huge volume spanning 0.0 < z < 1.2, and found only two 
[Oil] emitters with EW obs > 100 A and none with EW obs > 
120A. A higher fraction of EW obs > 120A [Oil] emitters 
have been found at higher redshift (e.g. Ajiki et al. 2006), 
but this is attributed to evolution rather than a selection 
bias at lower redshift. When this is taken into account, we 
expect to find ~0.21 ± 0.46 [Oil] emitters with such high 
EWs in our coverage of the FLS region. 

Once selected, each candidate must then be checked to 
confirm that it is a probable LAB; high proper motion ob- 
jects, passing satellites, objects in the haloes of saturated 
stars and amplifier crosstalk artefacts must all be eliminated. 
Thus, each candidate is inspected by eye to remove these 
"false positives" from the sample. The remaining candidates 
are considered robust LAB candidates and await follow-up 
spectroscopy. 



3 FIRST RESULTS - A LYMAN-a BLOB AT 
Z=2.83 

3.1 The Lyman-a emission 

Figure 1 shows our first spectroscopically confirmed LAB at 
17' l 09 m 52.3 s +59° 13' 21'.' 72 (J2000) in various wave-bands. 
It is clearly bright and extended in the Hell narrow-band 
filter, while it is comparatively faint and compact in the 
g' and i' band images. The Ly-a emitting structure is most 
clearly revealed in the continuum-subtracted figure 1(b) and 
has an observed equivalent width of EW obs = 657A. The 
photometric data for this new LAB at 2=2.83 are presented 
in table 1. We observed this LAB candidate with the William 
Herschel Telescope ISIS spectrograph with the blue and red 
arms using the R158B and R158R gratings respectively and 
with a 2" slit width. This results in a spectral resolution 
of ~ 15 A, and wavelength coverage between -3500-8300A. 
The one and two dimensional spectra can be inspected in 
figure 2, in which it is clear that the emission is extended 
over 12", corresponding to ~95 kpc at this redshift. 

We identify the line as Ly-a based on its EW, which 
is too high to be caused with any reasonable probability by 
a low-redshift interloper (see section [2j , particularly with 
the absence of any other emission lines in our wavelength 
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Figure 1. Observations taken in the various bands of the immediate vicinity of the LAB. The extended Ly-a emission is clearly visible 
in the continuum-subtracted and Hell narrow-band images shown in panes |l(b)| & |l(c)| Figure |l(c)| also has the slit positions for our 
follow up spectroscopy shown. Each pane is 26 X 26 arcsec 2 and North is up, with East to the left. The galaxy to the North and East of 
the LAB falls within the slit, and its continuum is apparent in figure [2(c) | It has been identified as being at zs!0.84 due to the presence 
of [Oil] 3727, [NeIV]2424, Mgll2799, and [CII]2326 emission lines in its spectrum. Figure fl (g) | has variable sky brightness across it since 
the LIRIS data have a chip boundary at the centre of the field of view; it is debatable whether or there is a detection in this band, and 
so has been included as an upper limit in table [TJ The u* band data are from Shim et al. (2006), and the R band data come from Fadda 
et al. (2004). 



range (3500 - 8300 A). The FWHM of the line is < 970 
km s" 1 , unresolved at our spectral resolution, and the to- 
tal Ly-a flux across the entire extent of the structure in 
the ISIS spectrum is 2.86 ± 0.07 x 10~ 16 erg s" 1 cm" 2 , the 
difference from the photometric measurement in the narrow- 
band filter being due to the geometry of the slit and the 4" 
aperture. Integrating across the entire extent of the LAB 
in the continuum-subtracted image, we estimate the total 
luminosity of the Ly-a emission to be 2.1 x 10 44 erg s _1 . 

The 2D spectra in figure 2 not only show a "clumpy" 
structure, but the emission line also seems to possess both 
a cut-off on the red side of the line (perhaps indicative of a 
rapid outflow velocity structure - see section f?] below), and 
hints of a possible shear which may be interpreted as rota- 
tion, although in neither case do we have the resolution to 
make a definitive statement; it is clear that higher resolution 
spectroscopy is required. 



3.2 The host galaxy 

LABs are reminiscent of the large haloes that surround pow- 
erful radio galaxies but with less than 1% of the associated 
radio emission (e.g. Reuland et al., 2003). Using the FLS- 
VLA survey at 1.4 GHz of Condon et al. (2003), we find 
that there is no detectable radio emission down to a radio 
flux-density limit of Si a ghz ~ 60/iJy (3<x). The LAB is also 
not detected in 610 MHz data taken using the Giant Metre- 
wave Radio Telescope (GMRT) from Garn et al. (2007) to 
a flux limit of SeiOMHz ~ 240/iJy (3<r). Thus it is clear 
that our LAB is significantly fainter than the radio lumi- 
nosity of powerful radio sources, with a radio luminosity of 
L1.4, ghz < 1.6 x 10 23 W Hz" 1 sr" 1 . 

At near-infrared wavelengths we obtained service-mode 
if-band data of the LAB using the UKIRT Fast Track Im- 
ager (UFTI; Roche et al., 2003) on source for 3780s [Fig- 
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Figure 2. One- and two-dimensional spectra of the LAB. The ID 
spectrum (top) was extracted over a region 4" in the direction 
SE to NW [see figure |l(b)| . For the 2D spectra (bottom left, 
and right), the spatial axis is horizontal and the spectral axis 
is vertical. Wavelength increases upwards and the images have 
been convolved with a Gaussian kernel with a standard deviation 
of 0.4". The Ly-a emission is extended over at least 12", and 
appears to have a potential cutoff on the red side in the left hand 
spectrum along with a slight shear, although the resolution is low. 
The right hand spectrum displays continuum emission to the left 
of the Ly-a line, which is due to the low-redshift galaxy to the 
North-East of the LAB, and not part of the LAB structure itself. 
The left hand spectrum was taken at position angle (PA) of 146° , 
while the right hand spectrum was observed at PA = 49° . 
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ure 1(h)]. From our observations we find that the LAB has 
Kab = 21.05lo'?5- We a l so obtained a 1200s J-band image 
from the Long-slit Intermediate Resolution Infrared Spec- 
trograph (LIRIS; Manchado et al. 1998) on the WHT, in 
which the LAB's continuum counterpart is not detected to 
a la magnitude limit of Jab > 19.74. 

The LAB falls just off the edge of the IRAC data in 
channels 1 & 3, but it is covered by channels 2 & 4 - right 
on the edge of the image. Although it looks by eye as though 
there is a faint detection in channel four, it is severely af- 
fected by reduced signal to noise near the chip edge and we 
only show the 4.5 /im data in figure 1. The 24 /j,m MIPS 
data however do cover the region around the LAB and a 
thumbnail image is displayed in figure l(j). There is a faint 
source within a 10" radius of the i'- and K-band contin- 
uum source, but it is considerably below the completeness 
limit of 300/iJy, and is firmly within the noise-dominated 
regime. It appears to occupy the same position as a knot of 
Lya emission to the South-East of the centre of the LAB 
in figure 1(c). If this 24^tm emission is real and associated 
with the knot of Ly-a in question, it would imply the pres- 
ence of a dust enshrouded burst of star formation. However 
deeper observations at this wavelength are clearly required 
to confirm this. It is considered unlikely that the association 
is indicative of an obscured AGN since it is offet from the 
continuum source in the centre of the LAB by ~ 50 kpc at 
this redshift, and there is also no other evidence of an AGN 
anywhere within the halo. 

Furthermore, we fit SED templates to the photometry 
of the embedded galaxy using our multi-wavelength data set 
(figure 3). We use the Bruzual & Chariot (2003) model SEDs 
normalised to the K-band flux to determine the mass of the 
embedded galaxy. The K-band flux is used for the normali- 
sation since it is the least likely to be affected by dust (the 
uncertainty on our 4.5/im point is prohibitively large), and 
is thus the most representative measure of the galaxy's old 
stellar emission available to us. We use a variety of different 
populations to fit synthetic SEDs to our observations; simple 
stellar populations of six different metallicities, twenty-one 
different ages between 10 Myr and 2.3 Gyr (the age of the 
Universe at 2=2.83 is 2.344 Gyr in our adopted cosmology), 
and initial mass functions (IMFs) from Chabrier (2003) and 
Salpeter (1955) are all used. Our LAB photometric data are 
poorly matched by a single simple stellar population (SSP), 
so linear combinations of young (< 200 Myr) and old (> 500 
Myr) SSPs of varying metallicities and IMFs are used in an 
attempt to better approximate the SED. 

We find that the best fit to our data points is a combi- 
nation of a young (20 Myr) and an old (1.0 Gyr) SSP with 
a Chabrier IMF and solar metallicity (Z = 0.02). In particu- 
lar, dust extinction is not required to give a good quality fit 
to the data, consistent with the absence of any significant 
24/im detection. With this SED fit, we calculate a stellar 
mass for this LAB of M ste uar = 2.79 ± 0.38 x 10 11 M , 
consistent with the most massive galaxies observed at high 
redshifts such as the hosts of powerful active galactic nuclei 
(e.g. Jarvis et al. 2001; McLure et al. 2006), and find that 
~20% of the LAB counterpart's stellar luminosity comes 
from the young starburst population, with the remaining 
~80% arising from the old population. 

We now consider the possibility that the counterpart 
galaxy within the LAB may be supplying the source of ion- 



ization for the LAB. We estimate from the model SED that 
the counterpart galaxy has a bolometric luminosity (blue- 
ward of the rest-frame Lyman limit, as would be expected 
for case B recombination) of ~2.4xl0 42 erg s _1 . Thus we 
conclude that at most only <1% of the LAB's Lyman a 
emission may be powered by the counterpart galaxy. 



4 DISCUSSION 

In this section we attempt to distinguish between the three 
ionization models for LABs discussed in section [T] 

In the absence of detections in both the 24^tm and 1.4 
GHz bands, the obscured QSO ionization model is highly 
unlikely. The flux limits in these bands essentially rule out 
the presence of a powerful AGN (both obscured and unob- 
scured). Moreover, the absence of any high-ionisation lines 
(such as CIV or CIII) indicative of an AGN, coupled with 
an optical/infrared SED (figure|3} that is approximated very 
well by two simple stellar populations, give no indication of 
the presence of an AGN at all. 

To estimate whether the central galaxy is a massive 
starburst galaxy we estimate its star-formation rate (SFR) 
from three different diagnostics. These are based on the 
ultra-violet, radio and Ly-a luminosities (see Brocklehurst 
1971, Kennicutt 1998, and Yun, Reddy & Condon 2001), 
where 

SFR(M yr' 1 ) = 1.4 x HT 28 L uv (erg s^Hz" 1 ) 

= 5.9 ± 1.8 x 10~ 29 Li.4GH*(erg s^Hz -1 ) 
= 9.1 x 10~ 43 L LyQ (erg s" 1 ). (3) 

For our LAB we find using the UV continuum luminosity 
(estimated from the SED fit) that the SFR of the embed- 
ded galaxy is 23.9 MQyr -1 . The radio data only provides a 
large upper limit (<12OOM0) so is not considered further. 
The SFR based on the Ly-a emission of the counterpart 
galaxy is more complicated, since it is difficult to separate 
the Ly-a emission of the galaxy from that of the halo. Thus, 
to estimate the SFR of the host galaxy we use the Lya lu- 
minosity measured in the central 4 arcsec which gives SFR 
~ 18 M©, consistent with our UV-based estimate. These 
optical SFR estimators assume the absence of dust; any ex- 
tinction would imply that these estimates represent a lower 
limit. However, we find no evidence for large dust masses 
using both the properties of our SED fit and the absence of 
any detectable reprocessed emission at 24/xrn. 

These results are highly significant since Fardal et 
al. (2001), using hydrodynamic simulations of the cold ac- 
cretion model, predict that LABs should have moderate 
star-forming galaxies at their centres. Although they pre- 
dict slightly larger SFRs of 60-600 Mq yr _1 , our values are 
within the right order of magnitude (particularly given the 
uncertainty in both the models and in the relations given in 
equations [3j . Furthermore, the cold accretion models from 
Dijkstra et al. (2006) also predict a sharp cut-off on the red 
side of the Ly-a line of the type hinted at by the red edge 
of our emission line in figure 3(a). 

On the other hand, the superwind model of Matsuda 
et al. (2004), estimates a SFR of 600 Meyr" 1 for a lOOkpc 
LAB, and draws a parallel between this and the submil- 
limeter detection (and thus high SFR ~ 1000 M Q yr _1 ) of 
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Figure 3. The SED of the galaxy immersed within the LAB 
(error bars), overlaid with the best fit redshifted simple stellar 
population models of Bruzual & Chariot (2003). The best fit con- 
sists of a combination of two stellar populations aged 30 Myr 
and 1.0 Gyr, both with metallicity Z = 0.02. ~20% of the LAB's 
bolomctric luminosity comes from the young population, with the 
remaining ~80% coming from the older population. 
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Figure 4. The continuum-subtracted surface brightness profile 
of the LAB from our narrow- and broad-band data. The flux- 
density is sky-subtracted emission line flux-density calculated in 
circular apertures, and the PSF of the image is shown with the 
dashed line. The solid line represents "model 2" from Dijkstra et 
al. (2006). 



LABI from S00 (Chapman et al., 2001). Our SFRs are con- 
siderably more quiescent; they may be able to cause some 
outflow, but it is probably unlikely that a starburst of this 
magnitude could cause superwinds of the types suggested by 
Taniguchi & Shioya (2000), Ohyama et al.(2003) or Matsuda 
et al.(2004). 

To further evaluate the likelihood of this LAB being a 
collapsing cloud of neutral hydrogen in a dark matter halo, 
we compare the Ly-a surface brightness profile (SBP) of 
the new LAB (figure 4) with that predicted by Dijkstra et 
al. (2006) . We find that there is a very good fit between our 
SBP and extended cold-accretion "Model 2" from Dijkstra's 
paper, with a total mass of 5.2 xlO 11 Mq. The biggest de- 
viation from Dijkstra's model is an apparent increase in the 
continuum-subtracted surface brightness around r ~ 50kpc. 
This is clearly not due to the nearby galaxy to the North- 
East of the blob centre as demonstrated in figure 1(b); it 
is rather a manifestation of clumpy sub-structure within 
the LAB. This is not predicted by the cold-accretion model 
of Dijkstra (2006), which assumes the gas within a LAB 



to be spherically symmetric. Future higher-resolution spec- 
troscopy over the whole halo would help to distinguish be- 
tween a starburst driven superwind and the cold-accretion 
model, however the current evidence favours the cold accre- 
tion model. 



5 CONCLUSIONS 

• We have presented our wide-field narrow-band survey, 
and demonstrated its ability to pick out LABs at high red- 
shift. Over the whole survey we expect to find a significant 
number of LABs to determine their nature fully. 

• We have shown the first confirmed LAB from our sur- 
vey as being extended over ~95kpc and having a Lyman a 
luminosity of 2.1 xlO 44 erg s _1 . 

• Our LAB is best described by the cold accretion ioniza- 
tion models of Fardal (2001) and Dijkstra (2006), although 
we cannot completely rule out the superwind model. It is 
highly unlikely that our LAB contains an AGN component. 
Further studies are required to strengthen our confidence in 
these assertions and also to gain better mass estimates for 
this LAB. 
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